Energy Loss of Leading Partons in a Thermal QCD Medium 
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We consider bremsstrahlung energy loss for hard partons traversing a quark-gluon plasma. Ac- 
counting correctly for the probabilistic nature of the energy loss, and making a leading-order accurate 
treatment of bremsstrahlung, we find that the suppression of the spectrum is nearly flat, with the 
most suppression at energies E ~ SOT (T the QGP temperature), in contrast to previous literature 
but in agreement with experimental data. This flat pattern should also be observed at the LHC. 
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Introduction 

In highly relativistic heavy ion collisions, production 
of hard partons precedes most other processes, sim- 
ply because the time scale of the production is short, 
T ~ 1/pt- In particular, the production of hard sec- 
ondary partons precedes the formation of a quark-gluon 
plasma (QGP). Therefore, the produced hard partons 
find themselves in an environment far different than the 
vacuum. The interaction between the parton and the 
environment influences the flnal spectrum of high 
hadrons in a non-trivial way. In particular, if the QGP 
is very dense, we expect energy loss, leading to an en- 
ergy dependent suppression of the high -pT spectrum, a 
phenomenon called "jet quenching." The extent of jet 
quenching can be used to learn about the QGP Ex- 
perimentally, the CERN SPS gave little evidence of jet 
quenching, but RHIC has seen a rather dramatic sup- 
pression of high -pT pions. 

The results from RHIC experiments^, 0,0,0 are sur- 
prising in many ways. At ~ 3 GeV, we already see 
a substantial suppression. Furthermore, the suppression 
continues more or less in constant fashion up to the high- 
est -pT measured so far (~ 10 GeV). The amount of sup- 
pression is also rather big; the ratio of high pT events 
to the number expected based on proton-proton data is 
1/5. 

Theoretically, it is well established that the main en- 
ergy loss mechanism of a fast parton is the bremsstrahl- 
ung of giuons in medium. The strength of the brems- 
strahlung in medium depends on a coherence effect called 
the Landau-Pomeranchuk-Migdal (LPM) effect. The 
first quantitative treatment of this effect (in QED) was 
by Migdal more recently it has been considered in 
QCD by several authors 0, 11,11113 ■ 

Jet quenching through bremsstrahlung energyloss has 
been considered by several authors before UM 0, 
IT^ IT^ H^ . The purpose of this paper is to revisit 
the energy loss calculation, particularly emphasizing two 
points. In many previous treatments, the (path length 
dependent) average energy loss is computed and applied 
to all primary partons. However, bremsstrahlung is not 
well described by an (path length dependent) energy loss. 
Bremsstrahlung energy loss is dominated by hard emis- 
sions; for any time scale up to the shower length, a parton 
typically has made about half its energy loss in a single 
emission. Therefore, if a sample of partons initially has 



the same energy, then after going through some path- 
length of medium, the distribution of final energies is 
as broad as the mean energy loss. This is illustrated in 
Fig. ^ This is especially important when the initial en- 
ergy distribution is a rapidly falling function of energy; 
the final distribution is dominated by the few partons 
which happen not to lose much energy. This has previ- 
ously been emphasized by Baier tt. al. '^ho found 
that it significantly influences the final spectrum (see also 
|l^). To account for it, we directly evolve the spectrum 
of the partons as they undergo bremsstrahlung energy 
loss. 

Second, most previous treatments have taken the LPM 
effect to be a parametrically large suppression. This is 
formally true whenever the parton and the emitted gluon 
are very energetic, E'parton, £^giuon > T. However, the 
LPM suppression (actual rate over Bethe-Heitler rate) is 
only a factor of 1/2 at E'giuon — lOT; for less energetic 
giuons the LPM suppression is small (and approxima- 
tions which take it to be large are badly in error) . When 
the parton spectrum is steeply falling, the most impor- 
tant bremsstrahlung events have Sgiuon ^ -Eparton; so 
at realistic parton energies one cannot take the LPM ef- 
fect to be parametrically large. Therefore, we will use the 
formalism developed by Arnold, Moore and Yaffe (AMY) 
flftj . which correctly treats the LPM effect (up to 0{gs) 
corrections) at all energies iJgiuon > 9sT. 

In this short paper we will concentrate on the qualita- 
tive features of the spectrum, specifically its shape (the 
Pt dependence of R) . The goal is to show that the trend 
seen in the data-that R falls slightly and levels off, but 
does not rise with px at least at accessible momenta-is 
the trend expected from bremsstrahlung energy loss. In 
previous literature it has appeared necessary to explain 
the data by invoking additional many body effects such 
as the Cronin effect and shadowing [l^, [13, uM ■ Par- 
ticular, we predict that even at LHC, the flat suppres- 
sion pattern seen in RHIC experiments should persist. 
Since we are only after these qualitative features, we will 
simplify the treatment somewhat, and consider a static 
thermal medium of quarks and giuons at a constant tem- 
perature T > Tc To be more quantitative, we need to 
take into account the nuclear geometry and the expand- 
ing system. We will do so in a subsequent publication. 



Brief Description of Formalism 
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We consider a small number of high energy partons 
traversing a thermalized QGP. The high energy partons 
are rare enough that their dominant interactions are with 
the thermal bath particles. We also work at leading non- 
trivial order in as- This is obviously an idealization, but 
it is hard to see how to do better at present. 

A parton traversing the QGP undergoes a series of soft 
scatterings with other constituents of the medium, with 
leading-order cross-section 



fsoft — Cs9s 



rC(qx). 



(1) 



(The group Casimir Cg is Cf=^ [quarks] or Ca—5 [glu- 
ons].) Here C{q±) is the differential rate to exchange 
transverse (to the parton) momentum q_L . In a hot ther- 
mal medium, its value at leading order in as is 
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qi(qi- 



■(27Ve+iVf). (2) 



These frequent soft scatterings can induce coUinear 
splitting (bremsstrahlung) of the parton. The time scale 
over which the parton and bremmed gluon overlap, in the 
absence of other scatterings, is r ~ ^ ~ , with x 

the momentum fraction of the gluon and pj_ the momen- 
tum component of the gluon perpendicular to the original 
parton. When (TsoftT is large, additional collisions typi- 
cally occur while the parton and gluon are still coherent; 
this can frustrate the original emission. 

This problem has been treated in the QCD context by 
BDMPS [3 and by Zakharov AMY have re-analyzed 
it with almost the same conclusions [Kij : we outline the 
physics and summarize their results. The probability of 
emission of a gluon of momentum k is schematically 



dk I 



(3) 



with the effects of soft collisions implicitly included in the 
time evolution. The time for the J • G insertion in the 
amplitude and its conjugate differ; to get the emission 
rate we must integrate over the difference of these times. 



dt' {p\ J^C^ {t')\p-k; k){p-k- kylGl (0) \p) . (4) 



dV 
dkdt 



The problem is then to evolve \p~k] k) {p\ between time 
and time t' . Its evolution equation is similar to a Boltz- 
mann equation, but with an extra phase accumulating 
term because the states \p) and |p— fc;fc) have different 
energies. Performing the time integration (taking the 
medium to be uniform on the scale of the formation 
time), the complete expression for the bremsstrahlung 
rate turns out to be. 



dT{p, k) 
dkdt 



Cs9l 



1 



1 



1 ± e-fc/^ 1 ± e-iP-^)/T 



a;3(l-a;)2 

xt 



3(l-a;)3 




2h-ReF(h,p, /c) , 



(5) 



where x= k/p is the momentum fraction in the gluon (or 
either quark, for the case g qq). The factors 1/(1 ± 
g-*:/^) are Bose stimulation or Pauli blocking factors for 
the final states, with — for bosons and + for fermions. 
h = p X k is the invariant describing the non-collinearity 
of the final states; h lives in a two dimensional transverse 
space. F(h,p, /c) is the solution of an integral equation 
describing how \p — k;k){p\ evolves with time, due to 
collisions and the energy difference of the two states; 



2h=idE{h,p,k)F{h) 



■C(qi) X 



(2^)2 

{{Cs ~ CA/2)[F(h) - F(h-fcqi)] 
+ (CA/2)[F(h)-F(h+pqi)] 
-|-(Ca/2) [F(h) - F(h- (p-fc) )] } ,(6) 



SE{h,p,k) = 



'■p-k 



2pk{p-k) 2k 2{p-k) 



p 
2p 



(7) 



Here m? are the medium induced thermal masses, equal 
to m|,/2 for a gluon and Cfg^T^/i = g^T'^/S for a quark. 
For the case of g qq, the (Cs — Ca/2) term is the one 
with F(h — pqx) rather than F(h — fcqj^). 

The treatment of BDMPS is the same, except that 
it uses ((j^-t-m^)^ in the denominator of Eq. ((SJ, and 
drops the mass terms in Eq. |(7J). These errors are not 
numerically significant. They also typically solve Eq. ® 
in a large h approximation, valid for large p/T, k/T but 
unreliable for k < lOT. 

Next, we use these expressions to evolve the hard gluon 
distribution Pg{p,t = 0) and the hard quark plus anti- 
quark distribution Pq{p,t = 0) with time, as they traverse 
the medium. The joint equations for Pq and Pg are 



dPqjp) 

dt 



dt 



PqiP+k) 



dT^^ip+k,k) dT^g{p,k) 
dkdt 



dkdt 



-2Pg{p+k) 



dTI^{p+k,k) 



dkdt 



Pq{P+k)- 



dT^g{p+k,p) 



dkdt 



-Pg{P+k)- 



dVSjp+Kk) 



dkdt 



-Pgip) 



drUp,k) , dT^Jp,k) 



dkdt 



dkdt 



-e{2k-p)] ,(8) 



where the k integrals run from — oo to oo. The integration 
range with fc < represents absorption of thermal gluons 
from the QGP; the range with k > p represents annihila- 
tion against an antiquark from the QGP, of energy (k—p). 
In writing Eq. ©, we used dT^g{p, k) = dT^g{p,p—k) and 
similarly for g — > qq; the function in the loss term for 



3 



9 ~^ 99 prevents double counting of final states. Since 
bremsstrahlung energy loss involves only small 0{gsT/p) 
changes to the directions of particles, these equations can 
be used for the momentum distributions in any particu- 
lar direction. For a single initial hard particle, they can 
be viewed as Fokker-Planck equations for the evolution 
of the probability distribution of the particle energy and 
of accompanying gluons. 

Numerical Evaluations 

We evaluate Eq. © by the impact parameter space 
method of [l^, on a grid of points p, fc; Eq. © is then 
solved on this momentum space grid by a second order 
Runge-Kutta algorithm. The errors are quadratic in the 
momentum discretization provided proper care is taken 
with the small k behavior of dT/dkdt. Discretization 
errors are under control and numerical costs are modest. 

To have a realistic pt spectrum for both peripheral and 
central collisions, we need a realistic initial px spectrum. 
Here we will use a parameterization taken from Ref. |l3j | , 



dN 



d^pj 



C 



(9) 



as the initial px spectrum of hard partons. Ref. |13| fitted 
this to ^/s = 200 GeV pp data, and obtained pq ~ 1.75 
GeV (consistent with the initial spectrum in Ref.|12j|). 

Results and Discussion 

Consider first the evolution of a sample of quarks, all of 
momentum p aXt — Q. Since bremsstrahlung is a statisti- 
cal process, the range of momenta broadens and the mean 
value falls. We show such an evolution in Fig. The unit 
time step in the figure was 16/f/^T, which is ~ 1/T for 
as = 1/3. The figure shows clearly that the momentum 
distribution broadens as fast as the mean value falls. For 
comparison, we also plotted the particle energy we would 
obtain if we took the energy loss process to be a steady 
dE/dt energy loss, determined by 



dE 

Energy loss approx.: — — — 

dt 



dk 



dkdt 



(10) 



in which case all the particles would have the same energy 
at any given time. We see that this is a bad description 
of the real energy distribution. 

In Fig|21 we show the ratio of the final and the initial 
quark spectrum calculated in two ways. The solid and 
dashed curves are the result of directly solving Eq.®, 
and the dotted and dash-dot curves are obtained by 
evolving particle energies according to Eq. (|10|) . As the 
quarks move with almost the speed of light, we need not 
distinguish the time and the length. With the plasma 
temperature of about 0.4 GeV, the energy range shown 
corresponds to about 6 GeV to 20 GeV, with a time range 
of 0.5 fm to 1.5 fm. 

From Fig|21 it is obvious that the dE/dt method is not 
a good approximation. When particles with a steeply 
falling energy spectrum lose energy, the final spectrum 
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FIG. 1: Time evolution of an initially monoenergetic ensemble 
of quarks. For comparison, vertical bars show the energy of 
the quarks if we take the energy loss to be a steady process, 
determined by Eq. IIOII . 
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FIG. 2: The ratios of the final and the initial momentum spec- 
tra for quarks. The black and red (solid) curves are calculated 
by solving Eq.JHJ directly. The blue and green (dashed) curves 
are calculated by first calculating dE/dt. 



is typically dominated by those few particles which hap- 
pened not to lose much energy. To give an extreme ex- 
ample, consider a spectrum falling as p~^^ , and com- 
pare the effect of two energy loss mechanisms. In the 
first, half the particles lose all their energy and half are 
unaffected; in the second, all particles lose exactly half 
their energy. In the former case, the spectrum is reduced 
by a factor of 1/2; in the latter case, it is reduced by 
(1/2)^" ~ 1/1000. Yet the mean energy loss of the two 
hypothetical processes is the same. The real case is some- 
where in between. Almost all particles lose energy, but 
a few particles happen to suffer fairly small energy loss, 
while others suffer large energy loss. This causes less sup- 
pression of the spectrum than if the energy loss process 
were uniform. 

This effect becomes more important as the spectrum 
of initial particle energies becomes steeper. We illus- 
trate this in Fig|31 which shows the ratio Pq{E,t = 
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FIG. 3: The ratios of the final and the initial momentum 
spectra for quarks. The solid curves are calculated by solving 
Eq.® directly. The dashed curves are calculated by first 
calculating dE/dt. The evolution time is 2/T. The integer n 
corresponds to having initial spectrum /o oc l/(po 



2/T)/Pq{E,t = 0) for initial distributions of the form 
C/{pI For small values of n, the dE/dt method 

is not that bad. However, as n grows, the result of the 
dE/dt calculation deviates more and more from the cor- 
rect one. 

Also note that the energy dependence of the suppres- 
sion factor is different from the dE/dt prediction; at 
p/T ^ 15 the suppression is actually becoming larger as 
p increases, though eventually this behavior turns over. 
This p dependence is similar to what is observed experi- 
mentally in the ratio plo t of AA and pp high px spectra. 
As noticed before [13, the previous LPM effect calcu- 
lations seem to be incompatible with the current RHIC 
data essentially because these calculations show rising ra- 
tios while the data shows a flat or slightly decreasing ra- 
tio. What we show here is that this discrepancy does not 
mean that the basic energy loss mechanism (LPM sup- 
pressed bremsstrahlung) is wrong. Instead, the failure 
of most previous approximations is largely due to mak- 
ing the dE/dt approximation. Part of the explanation 
also lies in the treatment of bremsstrahlung made here, 
in which the LPM effect is not taken as parametrically 



large, and absorption as well as radiation is allowed. 



Conclusion 

In this work, we have demonstrated that LPM sup- 
pressed bremsstrahlung can in fact explain the qualita- 
tive features of the high px experimental data. This is in 
contrast to previous findings. This is because of two fea- 
tures of our treatment. First, we determine the distribu- 
tion of final energies an energetic parton can end up with, 
rather than using the mean. This is important when the 
initial spectrum is steeply falling. Second, we do not as- 
sume that the LPM effect is parametrically large from the 
beginning, but use a treatment which handles the tran- 
sition between Bethe-Heitler and LPM correctly. Since 
this transition occurs at emitted gluon energies ~ lOT, 
such a treatment is necessary. We find that the ratio R 
of the data to the pp based expectations at first falls with 
energy, reaches a minimum around SOT, and then rises 
slowly thereafter. 

We have checked that R approaches 1 at very large px ■ 
However, this happens at px of a few hundred T . With 
a reasonable estimate of T = 0.5 GcV, this indicates that 
even at px — 25 GeV, the suppression pattern at LHC 
will be nearly fiat. This is in contrast to the previous 
finding of Gyulassy and Vitev [T^, who claim that the 
Cronin effect and nuclear shadowing are essential in de- 
scribing the fiat R from RHIC, but that these effects are 
small at the LHC, which will therefore observe rising R 
at large px- 

To turn our observations into a quantitative prediction 
for the Px spectrum will require a proper inclusion of the 
nuclear geometry and of the expansion and cooling of 
the QGP medium. We plan to address these issues in a 
future publication. 
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